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Long- Period Contributions to the Disturbing Functions of the 
Earth from the Seventh, Ninth, and Eleventh Zonal Harmonics 

by 

T.L. Felsentreger and W.J. Wickless 

Introduction w* y 

It is the purpose of this paper to present explicit formulas for the long- 
period terms due to the seventh, ninth, and eleventh zonal harmonics in the 
disturbing function of the earth in the case of an artificial earth satellite. The 
formulas are given for terms of the satellite’s orbital elements and the Delaunay 
variables. G. Giacaglia (1) has given general expressions for the long-period 
terms due to any of the zonal harmonics, which can be expressed in terms of 
the orbital elements. The apparent differences between the results of this paper 
and those in Giacaglia' s have been verified as due to the errors in the latter as 
it appears in the A.J. 


The contributions of the long-period terms to the mean motion of the argu- 
ment of perigee are also given. 


The Disturbing Function 

The earth's gravitational potential at a distance r from the center of the 
earth is 


u =JL 

r 


1 - tM X P. <■*> A 


I— n=2 
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where 


ft m GM 

G = gravitational constant 
M = mass of the earth 
R = radius of the earth 

J n = zonal harmonic coefficients (n = 2, 3, . . . ) 

P n = Legendre polynomials (n = 2, 3, . . . ) 

0 - geocentric latitude. 

Here, the earth’s radius is adopted as the unit of length. The seventh, ninth, 
and eleventh Legendre polynomials are 

P 7 (sin0) = -L. (429 sin 7 0-693 sin 5 0 + 315 sin 3 0-35 sin0) 

P g (sin0) =-L. (12155 sin 9 0- 25740 sin 7 0 + 18018 sin 5 0- 4620 sin 3 0 + 315 sin 0) 
1 2o 

P (sin0) =-i- (88179 sin u 0-23O945 sin 9 0 + 21879O sin 7 0-9OO9O sin 5 0 
11 256 

+ 15015 sin 3 0-693 sin0). 


Let 


a = semi-major axis of satellite’s orbit 
e = eccentricity of orbit 

i = inclination of orbital plane to equatorial plane 
i - mean anomaly 
f = true anomaly 
g = argument of perigee. 
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The Delaunay variables L, G, and H are 


L = -fjl a 
G = L Vl-e 2 
H = G cos i. 

Use is also made of the relations 


sin 4> = sin i sin ( f +g) 


a L * /I rx 

— = — (1 +e cos f). 
r G 2 


The long-period terms in the expansion of U as a Fourier series in l and 
g are given by 


1 
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Ud't, 


(see referenced) 


making use of the relation 


di =L il df. 


Denoting the long-period parts of U ? , U 9 , and U n by A ? F 2 A 9 F 2p , and 
A F , respectively, we have 

1 1 2p 


21/x 9 J ? e sini 
16384 L 3 G 13 


[lo(s-135— + 495 — -429— ) ( 33 -30— + 5 — jsin g 

|_ \ G 2 G 4 G 6 / \ L 2 L 4 / 

- 15(3-69— + 209 — - 143— ) (ll - 14 — + 3 — ^ sin 3g 

\ G 2 G 4 G 6 / V L 2 L 4 / 

+ 33 (l-l5— +27— -13— Vl-2— +— ^ sin 5 gl 
\ G 2 G 4 G 6 / \ L 2 L 4 / J 


3 


A n F, = - 

9 2p 
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A u F 2p =- 
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( u 2 u 4 116 118 U 10 \ 

1 -61 — + 570— -1802— +2261— -969 — ) x 
G 2 G 4 G 6 G 8 G 10 / 


X 


(■ 


323-680 £1 + 462 
L 2 


21-112«1 + 7SL) 

L 4 L 6 L 8 / 


sin 3g 


+ 7020 


^5 -265 — 

\ G 2 


+ 2130— - 5746 
G 4 


H 6 H 8 H 10 \ 

— + 6137— — 2261 — ) 

G 6 G 8 G 10 / 


X 


( 


G 2 G 4 G 6 

323-816 — + 678 — -200— +15 

L 2 L 4 L 6 


G^\ 

L 8 / 


sin 5g 


- 49725 


-4li! 1 250 "!-S142! + 437»!- 13351!} 

\ G 2 G 4 G 6 G 8 G 10 / 


X 


{ 


G 2 G 4 

19-60— +66 — -28 

L 2 L 4 


G 6 

L 6 



sin 7g 


+ 20995 /l -25 — + 90— -130— +85— -21— ^ x 
\ G 2 G 4 G 6 G 8 G 10 / 


L a G 2 , G 4 . G 6 G 8 \ . Q 

x [1-4 — +6 4 — + — sin9g 

V L 2 L 4 L 6 L 8 / 


Contributions to dg/dt 

Since the Delaunay set of variables is canonical with respect to the 
Hamiltonian F, which includes A 7 F 2p , A 9 F 2pl and A u F 2p , we have 

— = - — (see reference 2). 
dt 3G 
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Therefore, a computation of B(A .F 2p V3g(i - 7,9,11) provides the long-period 
terms in dg/dt due to the seventh, ninth, and eleventh zonal harmonics. 

The results are 
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